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The disorder parameter (/x) for the condensation of monopoles in QCD is estimated analytically 
in terms of gauge invariant field strength correlators. The continuum limit is discussed. 
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I. INTRODUCTION 

It is a long standing conjecture that confinement 
of color is due to dual superconductivity of QDC 
vacuumQQQ. 

Dual superconductivity means condensation in the vac- 
uum of magnetic charges, and can be detected by use of 
an order parameter (fi) , which is the vacuum expectation 
value of a magnetically charged operator fi 3 [3 [H . 
For gauge group SU (N) there are N — 1 different mag- 
netic charges and fi a ,(a = l,...,N — 1), are the corre- 
sponding /i operators. A gauge invariant definition of \x a 
can be given and the order parameters (/i a ) have been 
studied numerically on the lattice [8| [9j- The indepen- 

projectioi 

monopoles has been argued 0] E3 
numerically |l2j . 

Lattice simulations show that (/i a ) ^ in the confined 
phase and that (fi a ) = in the deconfined phase where 
magnetic charge is superselectedjLj. both for quenched 
@ |3 and for unquenched QCD The scaling be- 

havior of (fi a ) with respect to volume in the vicinity of 
the deconfining transition allows to determine the order 
and the critical indexes of the phase transition ^3 ■ 
In this paper I present an analytic computation of ()i a ) 
in terms of gauge invariant correlators of the gluon field 
strengths in the vacuum. The aim is to relate confine- 
ment to properties of the correlators, which have been 
extensively studied on the lattice^^] 0] 

The approximation that we shall use in our computa- 
tion is the same which is at the basis of the "stochastic 
vacuum" approach to QCD^^]0][2(|: we shall perform 
a cluster expansion of (p, a ) in terms of the correlators of 
the field strengths and we shall truncate it at the second 
cumulant, bilinear in the fields. We expect this approx- 
imation to work reasonably in studying infrared proper- 
ties , which involve correlations at large distances. 

In Sect. 2 we recall the definition and the properties 
of (/x a ) and we perform its cluster expansion. 

In sect. 3 we recollect what is known from lattice simu- 
lations on the two-point field correlators, and we compute 
the susceptibilities p a 



(1) 



which are convenient tools to describe dual superconduc- 
tivity and the deconfining transition 5] 6] . Some compu- 
tational details are presented in the appendix.. 

In Sect. 4 we discuss the results. Confinement and 
deconfinement are related to properties of the field cor- 
relators. We discuss in what respect the numerical deter- 
mination of the latter can be improved. As a byproduct 
we discuss the continuum limit of p a and of (fj, a ) and 
their renormalization. A similar behavior is known for 
the Polyakov line in quenched theory [2l| . 



II. CLUSTER EXPANSION OF (fi a ). 
A. The order parameter. 

In SU(N) gauge theories , with or without dynamical 
quarks, the operator which creates a monopole of mag- 
netic charge S- of species a ,(a = 1, ...,N— 1), at (x,t) 
is written as 000: 

p a (x,t)=exp[^- i J d 3 yb ± (x - y)Tr(^ a (y,t)E(y,t))} 

(2) 

g is the coupling constant, q is an integer. bj_(x — y) is 
the vector potential produced at y by a static monopole 
sitting at x, with 




H 



(3) 
(4) 



H is the magnetic field, including the Dirac string. 
Explicitely, putting the string along the direction n 



n A x 



x(x — n.x) 



(5) 



<& a (y,t) is a scalar field in the adjoint representation 
which has the form 221 



$ a (y,t) = U(y,tm iaq U\y,t) 



(6) 



with U a generic gauge transformation and ^ iag a diag- 
onal matrix of the form 



N - a N ~ a a 



a 



N N 
^(N-a) 



) (7) 
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U can be taken as a parallel transport from a ref- 
erence point xq to (jj, t) along a line C, U (x) — 



2 



Pexp[ig c A^dx^]. The gauge transformation U{x) 
identifies the abelian projection. 

In the gauge in which & a (y,t) is diagonal , (abelian 
projected gauge), the electric field which enters in eq(2) is 
the color component of the transverse electric field along 
the diagonal generator T a defined as 



T a =diag(0,0,. 



■1 , -1,0,0 0) (8) 

a , a + 1 



since Tr($ a ,T b ) = 5 ab . 

Ej_ is nothing but the conjugate momentum to Aj_ so 
that p a simply adds the field of a monopole of charge q 
to A a , . 



B. Cluster expansion. 



From eq(2) 
(**> = 



n! 



d 3 y v .. / d^nb^x-yi). 



Vl(x~y.M® a -Ey i (yi,t)...^ a -EYHyn,t)) (9) 

where we have put ($ a .E)(y, t) = Tr($ a .E(y, t)) in order 
to simplify the notation. 

The term in the expansion eq(9) linear in the electric 
field is zero by symmetry. In the "gaussian" approxi- 
mation of the stochastic vacuum model only two point 
correlators and the even terms survive. The combinato- 
rial factor is (2n — 1)!! so that 

(p a ) -exp[-^ J d 3 yi ..J dSyJiix-yx).. 

b^(x-y n )(($\Eyi(yi,tU$\Er(fo,t))] (10) 

or, changing variable to yi — x , using the translation 
invariance of the correlator, and the usual definition 



e= 2 4 

g 2 



we get 

(p a ) ~ exp[ 



162V 



(11) 



d 6 Vi I d^VKy^Kfc) 

K&Oh-ih)] (is) 



where <frf j is a gauge invariant correlator which coincides 
with the correlator of the field strength component along 
the color space direction a of eq(8) in the abelian pro- 
jected gauge. 

We shall identify <5"j (a?) with the gauge invariant cor- 
relator as measured on the lattice [lg]|I3- The lat- 
tice correlators are defined by connecting the two field 
strengths by a straight line parallel transport. In general 
the correlators do depend on the particular path used for 



parallel transport [23j . However in what follows we shall 
only use their qualitative behavior at short and at large 
distances, which should be independent of it. 

In the analysis of confinement it proves convenient to 
deal not with (p a ) itself, but with the related suscepti- 
bility defined by eq(l) . From eq(12) we get 



q 2 d 
WNd/3 



d 3 md 3 y2&i(yi)62(y 2 ) 

Kfrfo-fo)] ( 13 ) 



Eq(13) suggests a quantitative check of the gaussian 
approximation in the cluster expansion. Higher clusters 
would introduce terms proportional to higher powers of q 
. If p a proves to be proportional to q 2 to a good approx- 
imation higher clusters can be neglected. Old data on 
p a are consistent with that, but a careful and systematic 
numerical check is planned. 



C. p a from lattice. 



We recall that , since (p a ) = Z(S z + { g ) S) , (p a ) = 1 at 



(3 = and therefore 



(tx a ) = exp[ p{{3')dp 



(14) 



Lattice determinations of p a can be sum- 

marized as follows. If (3 C is the critical value of j3 cor- 
responding to the deconfinement transition, in the ther- 
modynamical limit , when the spatial volume L? s goes to 
infinity: 

1) For < f3 c p a tends to a finite limit , implying by 
use of eq(14) that (^ a ) ^ 0. 

2) For (3 > fa 



p ~ c — cL s 
c> 



(15) 



implying, again by use of eq(14) that (p a ) — in the 
thermodynamical limit. 

3) For (3 = (3 C p a obeys the scaling law HUH 



o x 



(16) 



with t = (1 — jr) the reduced temperature and v the 
critical exponent of the correlation length , A , of the 
order parameter. 



(17) 



We shall discuss how these features can emerge from 
eq(13). 
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III. COMPUTING p a 



and 



By general invariance arguments the electric field cor- 
relator at equal times has the form [ see e.g. eq(2.6) of 
ref. 13] 



- %[£(z 2 ) + \di{z 2 )\ + di^.D.iz 2 )] (18) 



with D,Di two invariant functions of k 2 . The depen- 
dence on a, (a = 1, ..N — 1) disappears in eq(12): this is 
in agreement with what is observed numerically^ 124| . 
Notice that this is true if we identify our gauge invari- 
ant correlators with the ones of the stochastic model as 
discussed above. A dependence on a , however, could in 
principle exist in our correlators. The last term in eq(18) 
does not contribute to the convolution with b± in eq(13) 
since = 0. After Fourier transform eq(13) becomes 

then 

;\fi / 7s4&i(*)6i(-*)*«(* a )] (19) 



with 



16 9/3 J (2tt) 



^ ij {e) = 5 ij [D{k 2 ) + \b l {k^] 



(20) 



Since k.b±(±k) = we can replace Sij in eq(20) by 

(21) 



k S -M^ihi and, putting 



1 ^ 



k 2 f(k 2 ) EE D(k 2 ) + -D^k 2 ) 



and making use of the equality 



(fc%- - kik^Wjj^ = H(k).H{-k) (22) 



with H (k) the Fourier transform of the magnetic field, 
we finally obtain 

P a = -f?4~^ I 4^\Hik)\ 2 f{k 2 )] (23) 



16 0/3^ J (2tt) 3 

If the Dirac string of the monopole is put along the di- 
rection n , say z, then 



b±(k) 



' A k 



k{k — k.n — ie) 



and 



H(k) = ~ik A b±(k) = - 



k 2 (ri.k) — ie 



\H(k)f 



1 1 



(24) 



(25) 



(26) 



q 2 d f d 3 k „,, 2 . , 1 1 sl . „. 

p= mwVJw m( e-it )] (27) 

At large (3 (high temperature) f(k 2 ) can be approximated 
by lowest order perturbation theory which is easily com- 
puted to be f(k 2 ) = , the dependence on (3 in eq(25) 
only comes from the explicit factor (3 and 



d 3 k 1 , 1 1 



16N J (2tt) 3 2k y k 2 k\ 



(28) 



The integral is easily computed [See Appendix A] with an 
ultraviolet cut-off - and an infrared cut-off j^— giving 



P = 



16N (2tt) 



■ [- V2L S + 2lnL s + constant] (29) 



Eq(29) means (/z) = as L s — > oo , or superselection 
of magnetic charge. This is indeed what is observed in 
lattice simulations @ Q 0] . Also the proportionality 
to q 2 of the slope with respect to L s is consistent with 
data0 . 

The field correlators D and D\ as measured on the 
lattice in the range of distances .1/m <-> 1/m are well 
fitted by a parametrization of the formal (2q|: 

D = Aoexp(~) + ^exp(~) (30) 

Aft X A a 

Di = A ieX p(~) + ^exp(~) (31) 

Af, X A a 

The slopes of both exponentials in D and D\ are equal 
within errors: Ah w -3/m, A„ w 2A&. 

The parametrization eqs(30)-(31) is inspired by the 
Operator Product Expansion (OPE). The terms propor- 
tional to ^4 correspond to the identity operator, the 
terms with only the exponential to the gluon condcn- 

2 

sate G2 : more precisely A + A\ = jgG2 with G2 — 

( Kb) fia pa \ 

\ g U [IV^* flU I 

In fact the parametrization eqs(30)-(31) cannot be cor- 
rect at short distances since the expansion of the second 
terms gives non-zero coefficients to terms in \ , \ , — im- 
plying the existence of condensates of dimension 1,2,3,, 
which do not exist. Similarly for the exponential terms 
the expansion of the exponential implies the existence of 
a condensate of dimension 5 , which again does not exist. 
The correct OPE should have the form 



D(Di) w^o -r + A + Cx 2 



(32) 



or better, to match the perturbative expansion at small 
distances 



1 



2 l (x + ie) 4 (x-ief 



}+A+Cx 2 + .. (33) 
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A possible simple form could be obtained by the replace- 
ments 

exp-(f ) — cxp-([(f ) 4 + e]* - e) (34) 



exp-(T-) — > exp -([(—; 

Afc At, 



(35) 



with suitable values of e, e' to match lattice data within 
errors. At large distances a strong infrared cut-off must 
exist in the confined phase , which should disappear at 
the deconfining transition. Such a cut-off implies for the 
term in A + 4^ that the space integral must be zero. A 
change of the form 



the singular part only comes from the first term in paren- 
thesis of eq(37) giving 



V2Kg 2 (3 
16A(2tt) 2 7 



(40) 



This is consistent with scaling eq(16). A quantitative 
numerical analysis is needed, however, to have a precise 
picture of the transition. 



Aexp(- 



^-m^-Yj (36) 



IV. CONCLUDING REMARKS 



is consistent with the lattice data on correlators within 
errors, and give zero integral. Parametrizations of the 
correlators which imply a change of sign at large distances 
exist in the literature [25|. 

Since the short distance behavior coincides with the 
perturbative one the singular term will give a result sim- 
ilar to that of eq(29), with L s replaced by the infrared 
cut-off A and a different value for the constant. The 
contribution of exponential term modified as in eq(34) is 
easily computed . The total result is 

Q 2 d , 1 , n k ,A, 

P = 77^77^ 7^2 -V^~ + 2M- + const.) 
loW op (27r) z a a 

--LjvA 6 3 A04+:U)] (37) 

From zero temperature up to say .95T C the correlators 
are constant, the only dependence on /3 is the explicit 
one in eq(27) so that 



I6A l (2tt 



1 , ^A , .A, 

(-V2— + 2ln(— ) + const) 



-—NXbA(A 
on 



(38) 



which is finite and volume independent, so that , by use 
of eq(14) (/i) ^ (dual superconductivity). 

This is true at fixed ultraviolet cut-off. In the contin- 
uum limit, however, when a — > p oc - and a renormal- 
ization is needed like for the Polyakov line |2lJ . Existing 
lattice data are consistent with this behavior ( see e.g. fig 
.2 of ref[8] ), but an extensive numerical check is planned. 

As T — > T c A — > co , the electric condensate (E 2 ) , and 
hence A and A\ vanish very rapidly |2(j [See e.g. Fig . 
10 of ref [26].] A a ,Ab do not change appreciably across 
T c . A diverges when approaching the critical point r = 0. 
Eq(37) must then be used. As a model example if 



A 

a 



Kln{r) 
-> 



(39) 



We have related the susceptibilities p a of the order pa- 
rameter for dual superconductivity to properties of the 
gauge invariant field strength correlators, by use of the 
stochastic vacuum model of QCD, i.e. by a cluster expan- 
sion of field strength correlators truncated at n = 2. A 
direct check of this approximation can be done by verify- 
ing that p a is proportional to the square of the magnetic 
charge. 

The validity of the perturbative approximation at large 
temperatures insures superselection of magnetic charge, 
i.e. that vacuum is normal. 

Dual superconductivity in the confined phase is related 
to an infrared cut-off of the two point correlator, which 
renders p a volume independent. 

The order parameter is anyhow ultraviolet divergent 
and needs a renormalization in the continuum limit, 
similarly to what happens to the Polyakov line in the 
quenched theory. 

The scaling behavior in the critical region depends on 
how the infrared cut-off diverges , and on how the electric 
condensate vanishes at the critical point. A direct mon- 
itoring of the large distance behavior of the two-point 
correlator is needed to confirm our assumptions, which 
go beyond the existing lattice data on correlators, as well 
as a more precise study of their behavior at the critical 
point. 
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APPENDIX A: A 

We compute the integral in eq(28) of the text. 

r= I (Ai) 



d 3 k 1 1 1 
(27r) 3 2fc ( fc2 " 



We integrate the first term in spherical coordinates, the 
second one in cylindrical coordinates. 



/ = 



dk 
~k 



dk. 



dk\ 



^2^ z+ kl) ] (A2) 
We have introduced as ultraviolet cut-off on each com- 
ponent of k the inverse of the lattice spacing a and as 
infrared cut-off the inverse of the spatial lattice edge L s a. 
Performing first the integral over k\ from k\min = a £ 2 

= l\ gives for the second integral in paren- 



thesis 



to A - " n 



k z min 



, , \ v* + kl)\ k A max (A3) 
-, k z min = -4—. Using the indefi- 

a 5 z aL s to 



Here again k z max 
nite integral 



dx 



c 2 ) = 



one finally gets 
1 



^ x2 + c2) +ln(x+^{x 2 +c 2 )) 
x 

(A4) 



(2tt) s 



[- V2L S + 2ln(L s ) + 2V3 -I - In 



(1 + Vz? 

2V2 



(A5) 



